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1. 1+3 SPACETIME DIMENSIONS AND THE METHOD OF SPHERICAL MEANS 

We would now like to derive an analog of d'Alembert's formula in the physically relevant case 
of 1 + 3 dimensions. As we will see, the analogous formula, known as Kirchhoff 's formula, can be 
derived through the following steps. 

• Given a solution u(t,x) to the 1 + 3 dimensional wave equation, we will define a spherical 
average of u centered at x. The average will depend on the averaging radius r. 

• For fixed x, we will show that a slight modification of the average will solve the 1 + 1 dimen- 
sional wave equation in the unknowns (t,r). With the help of our corollary to dAlembert's 
formula, we will be able to find an explicit formula for this modified function. 

• We will take a limit as the averaging goes to in order to recover an expression for u(t,x). 

This procedure is known as the method of spherical means. The final result will be stated and 
proved as a theorem. Before proving the theorem, we will develop some preliminary estimates. We 
will use spherical coordinates (r, 0,0) e [0,oo) x [0,7r) x [0,27r) on R 3 . Recall that if the spherical 
coordinates are centered at the Cartesian point {p l ) p 2 ) p 3 ) ) then the standard Cartesian coordinates 
(a; 1 ,^ 2 ,^ 3 ) are connected to spherical coordinates by 



(1.0.1a) x 1 = p 1 + r sin 9 cos 0, 

(1.0.1b) x 2 = p 2 + r sin 9 sin 0, 

(1.0.1c) x 3 =p 3 + rcos9. 

Also recall that the integration measure associated to -B r (0) is da = r 2 du, where du d = sin 9d9d<p. 
Here, oj represents the angular variables. We will abuse notation by using the symbol u> to denote 
both the angular coordinates (9, 0), and alternatively as the corresponding point (sin 9 cos 0, sin 9 sin 0, cos 9) e 
dB 1 (0). 

Proposition 1.0.1 (Spherical averages). Let u(t, x) e C 2 ([0, oo) x M 3 ) be a solution to the 1 + 3 
dimensional global Cauchy problem 



(1.0.2a) 
(1.0.2b) 
(1.0.2c) 



-d 2 u(t, x) + Au(t, x) = 0, (t, x) e [0, oo) 



u(0,x) = f(x), 
d t u(0,x) = g{x), 



X € 



X € 
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For each r > 0, define the spherically averaged quantities 



U(t,r; 


x) 


def 


1 

47rr 2 


F(r; 


x) 


def 


1 

47rr 2 


G(r; 


x) 


def 


1 

4 ?rr 2 



(1.0.3a) 
(1.0.3b) 
(1.0.3c) 

and their related modifications 

(1.0.4a) 
(1.0.4b) 
(1.0.4c) 



— - f u(t,a)da= — [ u(t,x + ru) dco, 

ir 2 JdBr(x) 47T JajeaBi(O) 

f f(cr)dcr, 

JdB r (x) 



,(x) 



U(t,r;x) d = rU(t,r;x), 



F(r;x) = rF(r;x), 



G(r\ x) = rG(r; x). 



ThenU(t,r;x) e C 2 ([0, oo)x[0, oo)) is a solution to the following initial + boundary -value problem 
for the one- dimensional wave equation: 



(1.0.5a) 
(1.0.5b) 
(1.0.5c) 
(1.0.5d) 

Furthermore, 



-d$U(t, r; x) + d 2 U(t, r; x) = 



(t,r)e[0,oo)x[0,oo), 
U(t,0;x) = 0, te[0,oo), 
U(0,r;x) = F(r;x), re(0,oo), 
d t U(0,r;x) = G(r; x), r e (0, oo). 



(1.0.6) 



\imU(t,r;x) = u(t,x). 

r^0 



Proof. Differentiating under the integral on the right-hand side of (1.0.3a), using the chain rule 
relation d r [u(t, x + ru)] duj - (Vw)(t, x + ruj) -udu = j^^u^t, a) da (where N(a) is the outward 
unit normal to dB r (x)), and applying the divergence theorem, we compute that 

(1.0.7) S r U - f V^uMio - JL_ f A,u(t,y)J>«. 



We now derive a version of the fundamental theorem of calculus that will be used in our analysis 
below. If h is a continuous function on M 3 , then using spherical coordinates (p,u) centered at the 
fixed point x, we have 

(1.0.8) 

9 r / h(y) d 3 y = d r / / P 2 h(p, x + pus) dudp = / r 2 h(r, x + rtu) dtu A = / h(a) da. 

JB r (x) ' J0 JweOBi(O) JtuedB^O) JdB r (x) 



Multiplying both sides of (1.0.7) by r 2 and applying (1.0.8), we have that 
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1.0.9) d r (r 2 d r U) = —d r [ \u(t,y) d 3 y = — [ Au(t,a)da 

AlT JB r (x) All JdB r (x) 



Differentiating under the integral in (1.0.3a) and using (1.0.2a), we have that 



1.0.10) &U(t,r;x) = — [ d?u(t,a) da = — — f Au(t,a)da. 

Anr 2 JdB r (x) Anr 2 JdB r (x) 



(x) 

Comparing (1.0.9) and (1.0.10), we see that 

1 



(1.0.11) 



dtU(t,r;x) = — d r (r 2 d r U) = d?.U(t } r;x) + -d r U(t,r;x). 



Multiplying both sides of ( l.O.llj ) by r and performing simple calculations, we see that 



;i.o.i2) 



dt[rU(t, r- x)] = d 2 [rU(t, r; x)]. 



(1.0.2b 


) - 


s ( 


1.0.5c) 



We have thus shown that the PDE (|1.0.5at) is verified by U d = rU. 



Using (1.0.2b) - (1.0.2c) and definitions (1.0.3b) - (1.0.3c), it is easy to check that the initial 



conditions (1.0.5c) - (1.0.5d) hold. Note that you will have to differentiate under the integral in 



(1.0.3a) in order to show that (1.0.5d) holds 



The limit (1.0.6) follows easily from the right-hand side of (1.0.3a), since u is continuous 



Finally, the boundary condition (1.0.5b) then follows easily from multiplying (1.0.6) by r before 
taking the limit r -> + . 

□ 

Corollary 1.0.2 (Representation formula for U(t,r;x)). Under the assumptions of Proposition 



1.0.1, for < r <t, we have that 



(1.0.13) U(t,r;x) = f rU(t,r;x) = -(f(t + t;x) - F(r -t;x 



p=r+t 



p=—r+t 



G(p; x) dp. 



Proof. (1.0.13) follows from (1.0.5a) - (1.0. 5d) and the Corollary to d'Alembert's formula. 



□ 
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